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ABSTRACT 

I have developed a procedure utilizing a Genetic-Algorithm-based optimization scheme to fit the ob- 
served light curves of an eclipsing binary star with a model produced by the Wilson-Devinney code. The 
principal advantages of this approach are the global search capability and the objectivity of the final re- 
sult. Although this method can be more efficient than some other comparably global search techniques, 
the computational requirements of the code are still considerable. I have applied this fitting procedure 
to my observations of the W UMa type eclipsing binary BH Cassiopeiae. An analysis of V-band CCD 
data obtained in 1994/95 from Steward Observatory and U- and B-band photoelectric data obtained in 
1996 from McDonald Observatory provided three complete light curves to constrain the fit. In addition, 
radial velocity curves obtained in 1997 from McDonald Observatory provided a direct measurement 
of the system mass ratio to restrict the search. The results of the GA-based fit are in excellent agree- 
ment with the final orbital solution obtained with the standard differential corrections procedure in the 
Wilson-Devinney code. 

Subject headings: binaries : close — binaries : eclipsing — binaries : spectroscopic 
methods:numerical — stars:individual (BH Cas) 



1. INTRODUCTION 

The problem of extracting useful information from a set 
of observational data often reduces to finding the set of 
parameters for some theoretical model which results in 
the closest match to the observations. If the constitutive 
physics of the model are both accurate and complete, then 
the values of the parameters for the 'best-fit' model can 
yield important insights into the nature of the object un- 
der investigation. 

When searching for the 'best-fit' set of parameters, the 
most fundamental consideration is: where to begin? Mod- 
els of all but the simplest physical systems are typically 
non-linear, so finding the least-squares fit to the data re- 
quires an initial guess for each parameter. Generally, 
some iterative procedure is used to improve upon this first 
guess in order to find the model with the absolute mini- 
mum residuals in the multi-dimensional parameter space. 

There are at least two potential problems with this stan- 
dard approach to model fitting. The initial set of parame- 
ters is typically determined by drawing upon the past ex- 
perience of the person who is fitting the model. This sub- 
jective method is particularly disturbing when combined 
with a local approach to iterative improvement. Many op- 
timization schemes, such as differential corrections (Proc- 
tor & Linnell 1972) or the simplex method (Kallrath & 
Linnell 1987), yield final results which depend to some 
extent on the initial guesses. The consequences of this sort 
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of behavior are not serious if the parameter space is well 
behaved — that is, if it contains a single, well defined mini- 
mum. If the parameter space contains many local minima, 
then it can be more difficult for the traditional approach to 
find the global minimum. 



2. GENETIC ALGORITHMS 

An optimization scheme based on a Genetic Algorithm 
(GA) offers an alternative to more traditional approaches. 
Restrictions on the range of the parameter space are im- 
posed only by observations and by the physics of the 
model. Although the parameter space so defined is often 
quite large, the GA provides a relatively efficient means 
of sampling globally while searching for the model which 
results in the absolute minimum variance when compared 
to the observational data. While it is difficult for GAs to 
find precise values for the set of 'best-fit' parameters, they 
are well suited to search for the region of parameter space 
that contains the global minimum. In this sense, the GA 
is an objective means of obtaining a good first guess for a 
more traditional method which can narrow in on the pre- 
cise values and uncertainties of the 'best-fit' . 

The underlying ideas for Genetic Algorithms were in- 
spired by Darwin's (1859) notion of biological evolution 
through natural selection. A comprehensive description of 
how to incorporate these ideas in a computational setting 
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was written by Goldberg (1989). In the first chapter of his 
book, Goldberg describes the implementation of a simple 
GA — involving several steps which are analogous to the 
process of biological evolution. 

The first step is to fill the parameter space uniformly 
with trial parameter-sets which consist of randomly cho- 
sen values for each parameter. The theoretical model is 
evaluated for each trial parameter-set, and the result is 
compared to the observational data and assigned a fit- 
ness which is inversely proportional to the root-mean- 
square residuals. The fitness of each trial parameter-set 
is mapped into a survival probability by normalizing to 
the highest fitness. A new generation of trial parameter- 
sets is then obtained by selecting from this population at 
random, weighted by the survival probabilities. 

Before any manipulation of the new generation of trial 
parameter-sets is possible, their characteristics must be 
encoded in some manner. The most straightforward way 
of encoding the parameter-sets is to convert the numeri- 
cal values of the parameters into a long string of numbers. 
This string is analogous to a chromosome, and each num- 
ber represents a gene. For example, a two parameter trial 
with numerical values x\ = 1 .234 and y\ = 5.678 would be 
encoded into a single string of numbers '12345678'. 

The next step is to pair up the encoded parameter-sets 
and modify them in order to explore new regions of pa- 
rameter space. Without this step, the final solution could 
ultimately be no better than the single best trial contained 
in the initial population. The two basic operations are 
crossover which emulates reproduction, and mutation. 

Suppose that the encoded trial parameter-set above is 
paired up with another trial having x 2 = 2.468 and y 2 = 
3.579 which encodes to the string '24683579'. The 
single-point crossover procedure chooses a random posi- 
tion between two numbers along the string, and swaps the 
two strings from that position to the end. So if the third 
position is chosen, the strings become 

123 145678 -» 123 1 83579 
246(83579 ^246 145678 

Although there is a high probability of crossover, this op- 
eration is not applied to all of the pairs. This helps keep fa- 
vorable characteristics from being eliminated or corrupted 
too hastily. To this same end, the rate of mutation is as- 
signed a relatively low probability. This operation allows 
for the spontaneous transformation of any particular posi- 
tion on the string into a new randomly chosen value. So if 
the mutation operation were applied to the sixth position 
of the second trial, the result might be 

24645|6|78 ^24645|0|78 

After these operations have been applied, the strings are 
decoded back into sets of numerical values for the param- 
eters. In this example, the new first string '12383579' 



becomes x\ = 1.238 and y\ = 3.579 and the new second 
string '24645078' becomes x 2 = 2.464 and y 2 = 5.078. 
Obviously, the new set of trial parameter-sets is related 
to the original set in a very non-linear way. This new 
generation replaces the old one, and the process begins 
again: the model is evaluated for each trial, fitnesses are 
assigned, and a new generation is constructed from the old 
and modified by the crossover and mutation operations. 
Eventually, after a modest number of generations, some 
region of parameter space remains populated with trial 
parameter-sets, while other regions are essentially empty. 
The robustness of the solution can be established by run- 
ning the GA several times with different random number 
sequences. 

Genetic Algorithms have been used a great deal for op- 
timization problems in other fields, but until recently they 
have not attracted much attention in astronomy. The ap- 
plication of GAs to problems of astronomical interest was 
promoted by Charbonneau (1995), who demonstrated the 
technique by fitting the rotation curves of galaxies, a mul- 
tiply periodic signal, and a magneto-hydrodynamic wind 
model. Many other applications of GAs to astronomical 
problems have appeared in the recent literature. Hakala 
(1995) optimized the accretion stream map of an eclipsing 
polar. Lang (1995) developed an optimum set of image 
selection criteria for detecting high-energy gamma rays. 
Kennelly et al. (1996) used radial velocity observations 
to identify the oscillation modes of a 5 Scuti star. Lazio 
(1997) searched pulsar timing signals for the signatures 
of planetary companions. Charbonneau et al. (1998) per- 
formed a helioseismic inversion to constrain solar core ro- 
tation. Most recently, Wahde (1998) used a GA to deter- 
mine the orbital parameters of interacting galaxies. The 
applicability of GAs to such a wide range of astronomical 
problems is a testament to their versatility. 

3. OBSERVATIONS 

I have applied a GA-based optimization scheme to my 
observations of the W UMa type eclipsing binary star 
BH Cassiopeiae. Due to an unfortunate historical ac- 
cident, this relatively bright object (my = 12.6) was ne- 
glected observationally for more than half a century prior 
to this investigation. 

3.1. Background 

The discovery observations of BH Cas were made by 
S. Beljawsky from the Simeis Observatory between June 
and September of 1928, and were later reported in the pa- 
per "34 New Variable Stars (Fourth Series)" (Beljawsky 
1931). The eleventh variable in his list was in the position 
of the star which later came to be known as BH Cas. Bel- 
jawsky gave it the temporary designation '353.1931 Cass' 
at the time of discovery. No maximum or minimum mag- 
nitudes were listed. 
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Follow-up observations of BH Cas were made by 
B. Kukarkin with the 180-mm Zeiss-Heyde refractor from 
the Moscow Observatory and/or the 1 80-mm Grubb-Mart 
refractor from the Tashkent Observatory. Kukarkin be- 
gan gathering observations in April 1936, and continued 
through the summer of 1937. The results of his observa- 
tions were reported in the paper "Provisional Results of 
the Investigation of 80 Variable Stars in Fields 6, 13, 15 
and 62." published in the journal Verdnderliche Sterne 
Forschungs-und Informationsbulletin (Kukarkin 1938). 
Based on his 62 visual observations and an examination of 
the plate collection at the Moscow Observatory, Kukarkin 
concluded that BH Cas was possibly W UMa type with a 
period near 0.5 days and an amplitude of 0.4 magnitudes. 

The most recent edition of the General Catalogue of 
Variable Stars (Kukarkin et al. 1969) contains a footnote 
on the entry for BH Cas which refers to a 1943 paper by 
P. Ahnert and C. Hoffmeister. After searching for BH Cas 
on photographic plates, and visually with a 350-mm re- 
flector, they concluded that "no star in the surroundings 
of the indicated place shows light variation..." (Ahnert & 
Hoffmeister 1943). After the appearance of this article, 
no observations of BH Cas appeared in the literature for 
more than 50 years. Considering the strength of their con- 
clusion, the absence of additional efforts to observe this 
object is not surprising. 

On the night of 19 February 1994, I obtained a series 
of CCD images of the region centered on BH Cas over 
a period of ~1.3 hours using the Spacewatch CCD at 
the Steward Observatory 0.91-m telescope on Kitt Peak. 
For each image, I measured the flux from BH Cas and 
from the comparison and check stars, GSC 01629 and 
GSC 01134 respectively. Plots of the relative flux over 
time revealed an increase in the brightness of BH Cas rel- 
ative to the comparison star while the latter remained con- 
stant relative to the check star (Metcalfe 1994). 

3.2. Photometry 

On twelve nights between September 1994 and October 
1995 I used the '2kBig' CCD at the Steward Observatory 
1.5-m telescope to obtain images of BH Cas and the sur- 
rounding area for photometric study. These observations 
allowed me to reconstruct a complete light curve in the 
V-band from a total of 432 data points, as well as partial 
curves in the R- and I-bands. 

For each night of observations, I reconstructed a flat- 
field image from a median combination of many data 
frames with non-overlapping star images. I used the 
IRAF ccdproc package to clean and calibrate each image, 
and the phot package to extract aperture photometry for 
BH Cas, the comparison star GSC 00784, and the anony- 
mous check star (see Figure 1). 

For ten nights in December 1996 1 used the 'P3Mudgee' 
3-channel photoelectric photometer (Kleinman, Nather & 
Phillips 1996) on the 2.1-meter telescope at McDon- 




FlG. 1. — A 10' x 10' image from the Digitized Sky 
Survey centered on BH Cassiopeiae. Comparison and 
check stars that were used in this study are circled: 
(1) GSC 01629, (2) GSC 01134, (3) GSC 00784, (4) 
anonymous check star for CCD data, and (5) GSC 00594. 

aid Observatory to obtain multi-color observations of 
BH Cas. These observations yielded complete U- and 
B-band light curves with 1041 and 1107 data points re- 
spectively. 

I followed the standard reduction procedure (Clemens 
1993), but some additional corrections were necessary. 
Intermittent problems with the filter wheel on all but the 
final night of observations required that I make small 
(~3%) normalization corrections derived from nightly 
phototube cross calibrations. The considerable color 
difference between BH Cas and the comparison star 
GSC 00594 (AB-V ~ 0.37) yielded significant second- 
order extinction for observations obtained at high airmass. 
I adopted reasonable values for the second order extinc- 
tion coefficients based on a study by Kim & Park (1993). 
In this paper, the authors demonstrated that the values 
of the extinction coefficients can change drastically from 
one night to the next, and that the 'second-order' coef- 
ficient could be anywhere from negligible to dominant. 
Based on their Table 2, 1 determined that the median value 
of the second-order coefficient in the B-band was 0.05. 
Adopting this value, the resulting corrections applied to 
the B-band light curve of BH Cas typically amounted to 
a few hundredths of a magnitude. Although the paper 
did not determine values of the second-order coefficient 
in the U-band, it demonstrated the large range of accept- 
able values relative to the first-order coefficient. This al- 
lowed me to adopt reasonable U-band coefficients based 
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TABLE 1 

Characteristics of Photometry Data Sets 



Telescope 


Julian Day 


Airmass 


Filter 


k> 


k" 


Steward 1.5-m 


2449767 


1.54-2.61 


V 


0.133 




Steward 1.5-m 


2449970 


1.12-1.26 


V 


0.079 




Steward 1.5-m 


2449971 


1.12-1.21 


v 


0.550 




Steward 1.5-m 


2449977 


1.12-1.43 


V 


0.362 




Steward 1.5-m 


2449978 


1.12-1.40 


V 


0.323 




Steward 1.5-m 


2449998 


1.12-1.19 


V 


0.147 




McDonald 2.1-m 


2450430 


1.14-2.81 


B 


0.253 


0.05 a 


McDonald 2.1-m 


2450431 


1.14-1.81 


B 


0.230 


0.05 a 


McDonald 2.1-m 


2450430 


1.14-2.81 


u 


1.254 


0.72 a 


McDonald 2.1-m 


2450431 


1.14-1.81 


u 


0.775 


0.33 a 



a Adopted value 



on the assumption that the corrected maxima should be 
approximately equal. The application of these corrections 
makes it impossible to determine reliable estimates of the 
O'Connell effect (Davidge & Milone 1984) from these 
data. Table [j] lists the airmass range and extinction co- 
efficients for all of the data sets used in this analysis. 

3.3. Spectroscopy 

On seven nights in September/October 1997 I obtained 
a time series of spectra of BH Cas using the Sandiford 
cassegrain echelle spectrograph (McCarthy et al. 1993) 
on the 2.1-m telescope at McDonald Observatory. I ad- 
justed the cross disperser and grating rotation to provide 
wavelength coverage from 5430 to 6670 A. The velocity 
resolution of this setup was ~2 km s -1 . 

I followed the standard IRAF reduction procedure for 
echelle spectra (Churchill 1995). Although most spec- 
tral orders contained at least one weak metal line, there 
were only two strong features: the Na D lines and Ha. 
The Na D lines were contaminated by narrow atmospheric 
and interstellar features, so the only line which offered the 
possibility of measuring reliable radial velocities was Ha. 

The velocities were derived from cross correlations of 
the observed spectra with a spectrum of the radial velocity 
standard star 31 Aql obtained on each night (y r = 100.5 
km s , Astronomical Almanac 1997). I determined the 
systematic stability of the instrument and reduction pro- 
cedure by cross correlating the spectra of 3 1 Aql obtained 
on different nights with each other. Systematic velocity 
shifts between the data sets ranged from 2-10 km s -1 with 
uncertainties in the range 1-3 km s . I corrected each of 
the derived velocities for these small nightly offsets. 

4. MODEL FITTING 

4.1. Wilson-Devinney Code 

Before about 1970, the standard approach to modeling 
close binary stars was that of Russell & Merill (1952). 



This geometrical model revolved a system of two similar 
ellipsoids to produce a light curve. While it was admit- 
tedly only a useful approximation of a true binary system, 
this model could be treated analytically, which was appro- 
priate for the tools available at the time. 

After a detailed treatment of close binary stars by Kopal 
(1959) which described the Roche equipotential surfaces, 
several authors attempted to calculate light curves by re- 
volving this physical model. Progress was limited by 
the computational facilities that were available at the 
time. Early attempts by Lucy (1968) and Hill & Hutch- 
ings (1970) represented substantial improvements over 
the Russell model, but they were still incomplete. Wil- 
son & Devinney (1971) introduced the model which has 
served as the foundation for many improvements over the 
past few decades (Wilson 1994). At present, the major- 
ity of published results for close binary stars are obtained 
using the Wilson-Devinney code (Milone 1993). 

There are two modes of operation in the Wilson- 
Devinney (W-D) code for overcontact binaries like 
W UMa systems. One of these modes uses a single contin- 
uous gravity darkening law for the entire common enve- 
lope to fix the temperature of the secondary star. The other 
mode keeps the secondary temperature a free parameter, 
allowing the model to be in physical contact without be- 
ing in thermal contact. Both overcontact modes require 
that the two stars have identical surface potentials, grav- 
ity darkening exponents, bolometric albedos, and limb 
darkening coefficients. Since the primary and secondary 
eclipses in the observed light curves of BH Cas have sig- 
nificantly different depths, I chose to use the overcontact 
mode with an adjustable secondary temperature. I used a 
blackbody radiation law and simple reflection. 

Some parameters of the model may be fixed initially 
from theory. I interpolated the values for the monochro- 
matic limb darkening coefficients from the tables in 
Al Naimiy (1978) to the proper temperature and wave- 
length. I used the effective wavelengths of the UBV 
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bandpasses and the temperature derived from the colors 
of BH Cas. The bolometric albedo (also known as the 
reflection coefficient) is usually set to 1.0 for radiative 
stars, and 0.5 for convective stars which transport some 
of the incident energy to other regions of the star before 
re-radiating it (Rucinski 1969). The dependence of the 
bolometric flux on local effective gravity is described by 
the gravity darkening exponent; a value of 1.0 corresponds 
to a direct proportionality between flux and gravity, and is 
appropriate for radiative stars (von Zeipel 1924). The flux 
at any point on the surface of a convective star is less de- 
pendent on the local gravity, and the exponent is thought 
to be 0.32 (Lucy 1967). Since the temperature of BH Cas 
is well within the convective regime, the gravity darken- 
ing exponents and reflection coefficients should be near 
the theoretical values for convective stars. These are the 
values that I assumed for all of the modeling. 

4.2. GAWD Code 

In the summer of 1997, I began to develop a simple 
GA-based optimization routine for the 1993 version of the 
W-D code (GAWD). Inspired by a plot in a paper by Stagg 
& Milone (1993) which showed a slice of ill-behaved pa- 
rameter space in the £l-q plane, I decided to try using a GA 
to optimize in this two-dimensional space first. The re- 
sults of this initial test were encouraging. I used the W-D 
code to generate a synthetic V-band light curve, and then 
I let GAWD try to find the original set of parameters from 
a uniform random initial sampling of 1000 points in the 
surrounding region of 0,-q space. After 10 generations, 
nearly 99% of the trial parameter-sets were statistically 
indistinguishable from the original set of parameters. 

The GAWD code quite naturally divided into two basic 
functions: using the W-D code to calculate light curves, 
and manipulating the results from each generation of trial 
parameter-sets. The majority of the computing time is 
spent calculating the geometry of the binary system for 
each set of model parameters. The GA is concerned only 
with collecting and organizing the results of many of these 
models, so I incorporated the message passing routines of 
the public-domain PVM software (Geist et al. 1994) to 
allow the execution of the code in parallel on a network of 
25 workstations. 

With this pool of computational resources, the GAWD 
code evolved to do more than I originally thought would 
be feasible. Adding the capability to fit for more than two 
parameters was simply a matter of extending the GA to 
deal with longer strings of numbers. The light curve mod- 
els take the same amount of time to run regardless of the 
number of parameters which are specified by the GA, so 
in some sense the extra parameters were added for free. 
This is not strictly true since increasing the dimensionality 
of the parameter space slows the convergence of the GA. 
After experimenting in the simple VL-q space for awhile, I 
added the capability to fit for the inclination i, the temper- 



ature ratio of the two stars T1/T2, and finally the temper- 
ature of the primary star 7\. Since the absolute tempera- 
ture of a star is unconstrained by observations in only one 
bandpass, I altered GAWD to fit light curves in the UBV 
bandpasses simultaneously. 

4.3. Global Search 

The first thing I had to do before starting GAWD was 
specify the ranges of the 5 parameters: (1) The measure- 
ment of the mass ratio (q = jj^) from the radial velocity 
data placed strong constraints on the allowed range. Ini- 
tially I fit a spectroscopic orbit to the radial velocity data 
allowing the orbital period to be a free parameter, and 
with the eccentricity fixed at 0.0. The period of the or- 
bit was consistent with the photometric period, so I fixed 
it for the final fit. I allowed GAWD to fit for a mass ra- 
tio between ±3cr of the final spectroscopic value. (2) The 
shape of the observed light curves indicates that BH Cas 
is overcontact, so I constrained the equipotential parame- 
ter SI to be somewhere between the inner and outer critical 
equipotentials for a given mass ratio. (3) The depth of the 
eclipses implies that the inclination is fairly high, but I al- 
lowed all values above i = 50°. (4) The temperature ratio 
is strongly constrained by the relative depths of the two 
minima, so I included all values between 0.93 and 0.97. 
(5) The photometric colors indicate that the components 
of BH Cas are fairly cool, so I allowed temperatures for 
the primary star between T\ = 4200 and 5000 K. 

After turning on the 25-host metacomputer with PVM, 
I started the GA master program on one of the faster com- 
puters. The process begins by reading the observational 
data into memory and assigning each light curve equal 
weight. A randomly distributed initial set of 1000 trial 
parameter-sets is generated, and each set of parameters is 
sent out to a slave host. Each slave calculates theoretical 
UBV light curves for the given set of parameters and re- 
turns the result to the master. Upon receiving a set of light 
curves, the master process sends a new job to the respond- 
ing slave and computes the variance of the three calculated 
light curves compared to the observed data. The fitness of 
the trial parameter-set is determined from a simple aver- 
age using the three individual variances. 

When the results of all 1000 trial parameter-sets are in 
hand, the master process normalizes the fitnesses to the 
maximum in that generation. One copy of the fittest trial 
parameter-set is passed to the next generation unaltered, 
and the rest of the new population is drawn from the old 
one at random, weighted by the fitness. Each trial is en- 
coded into a long string of numbers, which are then paired 
up for manipulation. The single point crossover operator 
is applied to 65% of the pairs, and 0.3% of the encoded 
numbers are altered by the mutation operator. The shuf- 
fled trial parameter-sets are decoded, and those which are 
still within the allowed ranges of parameters replace their 
predecessors. Computations of the new set of models are 
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distributed among the slave hosts, and the entire process is 
repeated until the fractional difference between the aver- 
age fitness and the best fitness in the population is smaller 
than 1%. As the evolution progresses, correlations be- 
tween parameters are revealed in the spatial distribution of 
trial parameter-sets after the worst have been eliminated, 
but before the final solution has converged. The parame- 
ters of the best solution in the final generation provide the 
initial guess where the traditional approach begins. 

At this point, it is instructive to ask: has the algorithm 
found the actual region of the global minimum? This is 
really an epistemological question. First, the algorithm 
performs a global random sampling of parameter space. 
After the initial convergence, the crossover and mutation 
operations continue to explore new regions of the param- 
eter space in an attempt to find better solutions. Repeating 
this procedure many times with different random number 
seeds helps to ensure that the minimum found is truly 
global; but no solution can, with absolute certainty, be 
proven to be the global minimum unless every point in 
the parameter space is explicitly evaluated. At best, an al- 
gorithm can only take steps to sample parameter space in 
a sufficiently comprehensive way so that the probability 
of converging to the global minimum is very high. 

The fitting routine supplied with the W-D code is a dif- 
ferential corrections procedure. This program calculates 
light and radial velocity curves based on the user-supplied 
first guesses for the model parameters, and then recom- 
mends small changes to each parameter based on the local 
shape of the parameter space. After the corrections have 
been applied, new curves are calculated and the shape of 
the local parameter space is again determined. New cor- 
rections are suggested, and the process continues until the 
suggested corrections for all parameters are smaller than 
the uncertainties. 

5. RESULTS 

I evolved the 'first guess' set of parameters for BH Cas 
using the GAWD code. The GA randomly populated the 
parameter space defined in §13 and allowed the set of 
trial parameter-sets to evolve. After 100 generations, the 
difference between the average set of parameters and the 
best set of parameters was insignificant. The parameter 
values in the final generation of 1000 trial parameter-sets 
averaged to (mean ± la): 



curve was assigned equal weight in the assessment of fit- 
ness. The radial velocity data were used only to restrict 



(i) _ 



[Ti/T 2 



0.474 ±0.002 
2.798 ±0.015 
69.°52±1.°42 
0.953 ±0.005 
4788 ± 106 K 



the range of possible mass ratios (see §|4.3[). 

For the final solution, I used the Wilson-Devinney dif- 
ferential corrections code with all of the radial velocity 
data, and a subset of ~100 data points from each light 
curve. Starting with the best solution from the GAWD code, 
I allowed Q, i, T 2 , and L\ to be free parameters. I itera- 
tively applied any significant corrections returned by the 
code to the parameters until all of the corrections were 
small relative to their uncertainties. Finally, I fixed T 2 and 
Q and allowed T\ and q to be free parameters instead. No 
significant corrections were returned by the code, and the 
final orbital solution for BH Cas yielded: 



n 

i 

T 2 



0.474 ±0.002 
2.801 ±0.003 
70?1±0.°2 
4790 ± 10 K 
4980 ± 10 K 



The uncertainties given here are probable errors from the 
DC output. In Figure 2, the best-fit model from the W-D 
code is shown with the data for comparison. The devi- 
ations of the fit from the data may be the result of the 
blackbody assumption for the radiation law, or due to 




Photometric Phase 



The GAWD code utilized only the light curve data to arrive 
at this result. All data points were included, and each light 



FlG. 2. — The best-fit model from the W-D code plotted 
with the observational data that were used for the fit. In 
the top panel, the observations are: U-band (upper), B- 
band (middle), and V-band (lower). In the bottom panel, 
the observations are: primary (solid points) and secondary 
(open points) components of BH Cas. 
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TABLE 2 

Properties of the Eclipsing Binary BH Cassiopeiae 



Binary Type 

Position a, 5 (J2000.0) 

Proper motion fi a , (arcsec century -1 ) . . . 

Period (days) 

Epoch of Primary Min. (JDq) 

Magnitudes at Max. my,mR,mi 

Spectral Type 

Distance (pc) 

Velocity amplitudes K\,K% (km s -1 ) 

Mass ratio q 

7-velocity (km s" 1 ) 

Space velocity (km s -1 ) 

Equipotentials Qi = 

Inclination i 

Temperatures T\ , (K) 

Masses m\,m2 (Mq) 

Pri. Radii r\ (pole, side, back) (Rq) 

Sec. Radii r2(pole, side, back) (Rq) 

Pri. Fractional Lum. Li(U,B, V) 

Sec. Fractional Lum. L2(U,B,V) 

Limb darkening x\ = X2 (U,B,V) 

Gravity darkening gi = g2 

Bolometric albedo r\ = r2 



W UMa (W-type) 
00 h 21 m 21: s 41, +59°09W2 
-0.8 ±1.0,-2.2 ±1.0 
0.405890(04 ±31) 
2449998.618(7 ±3) 
12.58, 12.32, 11.67 
~K4±2 
^160 

89.7 ±4.9, 189.0±4.9 
0.475 ±0.028 
-23.7 ±2.4 
^29.6 ±4.9 
2.801 ±0.015 
70?1±1.°4 

4790 ± 100, 4980 ±100 

0.74 ±0.06, 0.35 ±0.03 

0.98 ±0.01, 1.05 ±0.01, 1.12 ±0.01 

0.70±0.01,0.73 ±0.01,0.82 ±0.01 

0.414 ± 0.007, 0.399 ± 0.005, 0.385 ± 0.005 

0.586 ± 0.007, 0.601 ± 0.006, 0.615 ± 0.005 

1.000, 0.976, 0.805 

0.32 

0.50 



the assumptions made during the reduction procedure, or 
both. Copies of the data shown in Figure 2 are available 
in the electronic edition of this paper. 

I have listed the properties of BH Cas (with \a errors) 
in Table ^[ The position and proper motion were obtained 
from a pair of 51 -cm Astrograph plates on the measur- 
ing machine at Lick Observatory (A. Klemola, personal 
communication). The results were noted to be somewhat 
worse than is normally expected from Astrograph plates 
due to the weakness of the exposures and the location of 
the image near the plate edge. 

From 12 times of minimum light spanning nearly 2000 
cycles (Metcalfe 1997), I derived the following ephemeris 
for BH Cas: 

Mini = JD 2449998.618(7 ±3) + 
+ 0?405890(04±13)xE 

The residuals of this fit exhibit no systematic trend, and 
the single largest departure is ~ 1 x 10~ 3 days. 

I derived the magnitudes of BH Cas at maximum light 
in the V-, R-, and I-bands from CCD observations us- 
ing the Selected Areas for zero-point calibration (Landolt 
1992). The colors of W UMa stars change only slightly 
with orbital phase because the common envelope forces 
the secondary star to be hotter than it would be if iso- 
lated. Observationally, this appears to be true regardless 



of the mass ratio of the system (Shu & Lubow 1981). The 
V — I and R — I colors imply an effective temperature for 
BH Cas near T e ff = 4600 ±400 K, corresponding to a Main 
Sequence spectral type of roughly K4±2 (Bessell 1979). 

An isolated Main Sequence star with this same effec- 
tive temperature would have an absolute visual magnitude 
My = +7.0 ±0.7 (Allen 1973). When this is combined 
with the observed magnitude it provides an upper limit on 
the distance. The primary (more massive, cooler) com- 
ponent of BH Cas has a structure which is approximately 
the same as an isolated Main Sequence star of comparable 
mass. The same cannot be said of the secondary (less mas- 
sive, hotter) component. As a consequence, it is best to 
use the observed magnitude at the orbital phase when es- 
sentially all of the light is coming from the primary com- 
ponent. For BH Cas, this occurs during the minimum of 
the deeper eclipse when my = 13.1. The resulting upper 
limit on the distance is d ^ 160 pc. 

The mass ratio and center of mass radial velocity of 
the binary system listed in Table || comes directly from 
the spectroscopic orbital solution. I obtained identical re- 
sults whether or not the few observations during eclipses 
were included in the fit. The transverse component of the 
space velocity is determined by combining the distance 
with the measured proper motion. The total proper mo- 
tion is /i tot = (^q+Z^) 1 / 2 = 2.34± 1.0 arcsec century" 1 . At 
a distance of 160 pc, this motion corresponds to a trans- 
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verse velocity of vj - 11 .1 ±7.6 km s . Combining this 
with the measured radial velocity of the system, the total 
space velocity is vs <J 29.6 ±4.9 km s" 1 . Thus, BH Cas 
appears to be a member of the disk population. 

In addition to the mass ratio, the spectroscopic orbit also 
yields the total mass of the binary system in terms of the 
semi-major axis, the orbital period, and some fundamen- 
tal constants. If the semi-major axis is expressed in terms 
of the mass ratio, K amplitudes, and inclination, then the 
individual masses can be calculated. I list the masses ob- 
tained by using the best-fit value of the inclination. 

Once the mass ratio is established and the scale of the 
system is determined from the K amplitudes, the best-fit 
value for the equipotential parameter leads directly to the 
absolute radii of the two stars. Since BH Cas is an over- 
contact system, the stars are tidally distorted, and no sin- 
gle number can adequately define the radii; I list the radii 
at the pole, side, and back. 

6. DISCUSSION 

The orbital period of BH Cas places it comfortably 
within the range of periods typical of many W-type 
W UMa systems. The mass ratio is near the average for 
this class, but the individual masses are near the low end 
of typical values. The mass of the secondary star seems 
particularly low, but this is not uncommon among W UMa 
systems. The secondary mass in V677 Cen, for example, 
is 0.15 Mq (Barone, Di Fiore & Russo 1993). In a recent 
paper containing a large self-consistent sample of abso- 
lute elements for W UMa systems (Maceroni & van't Veer 
1996), fully 12% of the W-type, and 18% of all W UMa 
secondary stars had masses less than or equal to the sec- 
ondary mass derived for BH Cas. The total system mass, 
however, is smaller than that of any W UMa star in this 
same sample, the smallest of which is CC Com (1.20 Mq 
total). 

The absolute radii of the components of BH Cas are 
larger than expected for stars with such low masses. The 
average radii are similar to those of TY Boo, another W- 
type system which has roughly the same mass ratio as 
BH Cas, but a systematically higher total mass (Milone et 
al. 1991). There are a number of possibilities which could 
have led to biased estimates of the masses of BH Cas. I 
derived the velocities exclusively from the Ha line since it 



was the only strong feature available. I used an unbroad- 
ened radial velocity standard for the cross correlation tem- 
plate. Also, I assumed Gaussian profiles for the correla- 
tion peaks. 

The fit to the radial velocity data could also be a prob- 
lem. There is considerable scatter in the observations, 
perhaps enough to accommodate the increased K ampli- 
tudes that would be needed to make BH Cas as massive 
as TY Boo. A set of high dispersion spectra of BH Cas 
obtained at the Multiple Mirror Telescope in 1994 has re- 
cently been re-analyzed using two dimensional cross cor- 
relation techniques with synthetic templates (G. Torres, 
personal communication). The re-analysis yielded a mass 
ratio consistent with the value derived from the echelle 
spectra, but the best-fit K amplitudes were considerably 
larger. The correlations which led to these results, how- 
ever, were very weak. 

A significant x-ray signal at the position of BH Cas was 
serendipitously discovered by Brandt et al. (1997). They 
found an x-ray flux of 4.2 x 10~ 14 erg cm -2 s" 1 . The im- 
plied upper limit on the x-ray luminosity from the distance 
given above is L x <J 1.0 x 10 28 erg cm -2 s -1 . The detection 
of x-rays from BH Cas is not too surprising considering 
that most W UMa systems show clear signs of magnetic 
activity (Maceroni & van't Veer 1996). The mechanism 
for x-ray production is generally thought to be related 
to magnetic reconnections in the corona (Priest, Parnell 
& Martin 1994). The value of the x-ray luminosity of 
BH Cas is lower by 1-2 orders of magnitude compared 
to other nearby W UMa stars (McGale, Pye & Hodgkin 
1996). 
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nificantly improved the manuscript. This work was sup- 
ported in part by a fellowship from the Barry M. Goldwa- 
ter Scholarship and Excellence in Education Foundation, 
an undergraduate research grant from the Honors Cen- 
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